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1 .  INTRODUCTION 


In  an  earlier  examination  of  the  Cerenkov  radiation  of  charged 
particles  passing  over  a  dielectric  plate,*  it  was  noted  that,  under 
certain  conditions,  a  beam  of  these  particles  could  support  a  number  of 
intrinsic  electromagnetic  modes.  This  was  a  surprising  result,  as  we 
had  not  before  encountered  a  discussion  of  these  modes,  and  a  search  of 
the  available  literature  turned  up  only  one  study  of  these  modes,  with  a 
brief  analysis,  for  slab  geometries  only,  of  the  results.1  In  view  of 
the  relevance  of  this  topic  and  our  findings  to  the  expanding  field  of 
near-millimeter  waves,  we  decided  to  continue  our  analysis. 


This  report  discusses  the  derivation  of  the  interaction  between  the 
electromagnetic  field  and  the  electron  beam.  Although  this  derivation 
may  be  found  in  many  books  containing  sections  on  Cerenkov  radiation  or 
monographs  on  travelling-wave  devices,2'3  the  derivation  appears  here  to 
familiarize  the  reader  with  our  particular  notation  and  assumptions. 
This  derivation  generates  the  wave  equations  satisfied  by  components  of 
the  electromagnetic  fields  within  a  suitably  restricted  electron  beam. 
The  theory  applies  to  the  case  of  a  ribbon  type  of  electron  beam  in  free 
apace.  It  is  demonstrated  that  the  form  of  the  wave  equation  for  this 
situation  suggests  a  superficial  resemblance  between  the  propagation 
constants  of  an  electron  beam  and  those  of  a  dielectric  of  identical 
geometry.  By  selecting  some  realistic  parameters  (such  as  current 
density  and  beam  size),  we  find  a  few  modes  of  ribbon  beam. 


1Xstvan  Palocz,  A  Leaky  Wave  Approach  to  Cerenkov  and  Smith-Purcell 
Radiation,  Ph.D.  Dissertation,  Polytechnic  Institute  of  Brooklyn, 
University  Microfilms,  Inc,,  Ann  Arbor,  MI,  # 62-5658  (1962). 

2M.  Chodorow  and  C.  Susskind,  Fundamentals  of  Microwave  Electronics , 
McGraw-Hill  Book  Co.,  New  York  (1964),  141-177. 

3R.  G.  E.  H utter.  Beam  and  Wave  Electronics  in  Microwave  Tubes,  D.  Van 
Nostrand  Co.,  New  York  (1960),  182-233. 

*Clyde  A.  Morrison  and  Richard  P.  Leavitt,  Cerenkov  Radiation  from  an 
Electron  Beam  Passing  over  a  Dielectric  Slab  Backed  by  a  Metal  Surface, 
Harry  Diamond  Laboratories  (draft). 


This  same  procedure  is  then  used  to  find  the  modes  of  a  cylindrical 
beam  in  free  space.  Because  no  available  literature  deals  with  this 
case,  we  examine  these  modes  in  considerable  detail  and  discuss  the 
amount  of  energy  coupled  out  of  these  modes. 


2.  ELECTRON  BEAM  WAVE  EQUATIONS 

These  are  Maxwell's  equations  listed  in  the  form  used  in  this 
report : 


V  x  h  =  -ikE  +  (4ir/c)Jt  ,  (1) 

V  x  E  =  ikH  ,  (2) 

V  •  H  =  0  (3) 

V  •  E  »  4ltPt  ,  (4) 


*►  +  —  i  t" 

where  H  and  E  vary  as  e  and  k  =  u)/c .  In  equations  ( 1 )  to 

+  +  > 

(4),  H  and  E  are  the  magnetic  and  electric  fields,  is  the  current 

density,  and  is  the  charge  density.  Unrationalized  cgs  units  are 

used  throughout  this  report.  We  make  the  following  linearizing 

approximations : 

Pt  “  Po  +  P  '  (5) 


6 


(6) 


vt  =  VQez  +  V  , 


Jt  -  Joe2  +  J  • 


*  + 

where  pQ,  VD,  and  JQ  are  constant  in  time  and  space  and  p,  V,  and  J 
as  e-i<*lt.  We  assume  that  p,  V,  and  J  are  small  in  comparison  with 
VQ,  and  JQ.  Also, 


=  Ptvt 


The  Lorentz  force  equation  is  used  here  in  the  form 


Here, 


and 


Y  -  1/(1  -  v|/c2)l/2 


d_  _  9 _ 

dt  =  at 


+  V. 


$  . 


Prom  equations  ( 1 )  and  ( 4 )  comes  the  conservation  condition 

p  +  $  •  J  =  0 


(7) 


vary 

“o' 


(8) 


(9) 


(10) 
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Inserting  equations  (5),  (6),  and  (7)  into  equation  (8)  and  keeping  only 
first  order  terms,  we  obtain 


^o  *  Povo®z 


(ID 


J  *  pe2VQ  +  p0V 


(12) 


Note  that  all  variables  in  equation  (11)  are  constant  with  respect  to 
time  and  position. 

The  next  step  is  to  search  for  travelling-wave  solutions  of  the 
fundamental  equations.  We  assume  that,  for  all  functions  of  z, 

.  iaz 

F(x,y,z)  =  F(x,y)e 


Using  equation  (9)  and  the  above  assumptions,  the  following  expressions 
are  found: 


V 

x 


V 

y 


v 

z 


ief  E  -  8H  ) 
mY3(u  -  oV) 

ie(E  +  6H  ) 
my3  (u)  -  aV) 

ieE 

_ z 

nry3(u)  -  oV) 


(•3) 


(14) 


(15) 


where  8  *  V/c. 
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By  substituting  equations  (13),  (14),  and  (15)  into  equation  (12), 
we  get  these  results  for  the  components  of  J : 


iui2f  E  -  6H  ) 
x  y' 


j  =  — * — - - * 

X  4tty3(u)  -  aV) 


i(D2  f  E  +  BH  1 

j  =  -  p  y . JL 

y  4hy3  (a)  -  aV) 


4tty3(ii)  -  aV) 


where  ui2  =  -4im_e2/m.  To  derive  equation  (18),  it  is  necessary  to  use 
P  ° 

equations  (16,*  and  (17)  in  equations  (10)  and  (4).  Doing  so  gives  also 


igui2E 

_ E_5 _  . 

4ityc(w  -  aV) 


We  can  now  get  the  wave  equation  for  our  beam  by  taking  the  curl  pf 
equation  (1)  and  using  equations  (2),  (16),  (17),  and  (18): 


V?H  +  k2  -  a2  - 


H  =  0 


and  similarly  for  E, 


72E  +  ^1 


k2  -  a2  - 


where 


Except  for  a  change  in  propagation,  these  equations  are  the  same  as 
one  would  find  for  an  ordinary  medium.  There  are  two  sets  of  solutions 
corresponding  to  transverse  magnetic  (TM)  and  transverse  electric  (TE) 
modes . 


3.  RESTRICTED  BEAM 

Consider  now  an  electron  beam  restricted  such  that  V  =  0,  V„  =  0, 

x  y 

and  as  before  (eq  15)  we  have 


V 

z 


-ieE 

z 

=  -■  —  — 
my3(u)  -  oV) 


(22) 


The  assumption  that  Vx  =  Vy  =  0  is  physically  realizable  with  the 

application  of  a  strong  magnetic  field  in  the  z-direction. 


Also,  Jx  *  0,  Jy  =  0,  and  from  continuity  (wp  =  aJz)  we  obtain 


J 

z 


iu)2uE 

_ P  « 

4by3{u  -  aV)2 


(23) 
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and 


P 


iai»)2E 

_ E_5 _ 

4irr3  (ui  -  aV)2 


(24) 


The  restriction  on  the  beam  has  a  considerable  effect,  as  can  be  noted 
by  comparing  equation  (23)  with  equation  (18)  and  equation  (24)  with 
equation  (19).  Using  equations  (1)  to  (4),  we  can  find  the  wave 
equation 


$($  •  E)  -  V2E  =  k2E  +  ik  —  J 

c 


(25) 


which  can  be  written  in  component  form  as 


|-  (4irp)  -  V2E  =  k2E 
3x  xx 


(26) 


f-  ( 4irp )  -  V2E  = 
3y  y 


k2E 


(27) 


I-  (4xp)  -  V2E  *  k2E  +  ik  —  J 
3z  z  z  c  z 


(28) 


Substituting  equations  (23)  and  (24)  into  equation  (28)  yields 


(29) 
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Vi 


where  kz  =  idz/cz.  This  equation  is  frequently  encountered  in 
P  P 

discussions  of  travelling-wave  devices.2  The  factor  k2/y2(k  -  aB)2 

gives  rise  to  solutions  for  a  near  nB  =  k,  in  which  the  phase  velocity 
of  the  wave  is  nearly  equal  to  the  velocity  of  the  electrons  in  the 
beam.  The  beam  must  be  restricted  in  order  to  generate  these  additional 
solutions.  No  such  behavior  is  possible  from  the  wave  equations  for  the 
unrestricted  beam  given  in  equations  (20)  and  (21). 

3.1  Ribbon  Beam 


where  the  y-axis  is  perpendicular  to  the  paper.  Assume  the  beam  to  be 
infinite  in  the  y-  and  z-directions  and  3/9y  =  0.  Then,  in  component 
form.  Maxwell's  equations  become 

-aE  =  ikH  ,  -ioH  =  -ikE  ,  (31a,b) 

y  x  y  x 

3E 

ioEx  =  ikHy  +  ~  ,  iaHx  =  -ikEy  ,  (32a,b) 

2M.  Chodorow  and  C.  Susskind,  Fundamentals  of  Microwave  Electronics , 
McGraw-Hill  Book  Co.,  New  York  (1964),  141-177. 


(33a, b) 


3_ 

3x 


E 

y 


“ikHz 


|-H 

9x  y 


*  ikE 

z 


+  iaE 

z 


4rtp 


¥x  Hx 


iaH 

z 


(34a, b) 


Further,  considering  TM  modes  only,  Hz  =  0,  which  from 

equations  (31)  to  (34)  immediately  gives  E  =  0  and  H  =  0.  From 

y  « 

equations  (31b)  and  (32a), 


for  all  values  of  x. 


H 

y 


ik 


k2  -  a2 


3_ 

3x 


E 

z 


Inside  the  beam,  equation  (30)  becomes 


-  E  +  (k2  -  a2)  1 

3x2  2 


k2  1 

P  . 

Y3 (k  -  a&)2_ 


E 


z 


0 


and  outside  the  beam  =  0  so  that 


32  9  ? 

-  E  +  (k2  -  a2 )E  =  0 

3x2  Z 


(35) 


(36) 


(37) 


The  symmetry  of  this  geometry  allows  one  to  choose  either  even 
or  odd  solutions  to  these  equations;  thus,  we  can  restrict  ourselves  to 
x  >  0.  Looking  at  the  even  solutions,  then,  we  can  write 


Ez  =  A  cos  ^k^x^e^aZ 


0  <  x  <  b/2 


(38) 
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VI 


^fi**-*f>*'  ■' 


1 


( 


and 


-k^  A  sin 


(¥) 


ik^Be 


ik2b/2 


(45) 


Dividing  equation  (45)  by  equation  (44)  gives  the  following 
transcendental  equation : 


-k 


1 


(46) 


This  expression  has  no  solutions  if  both  k1  and  k2  are  real. 
Assume  that  k,  is  real  and  k2  =  ik',  where  k£  is  real.  Thus,  k^  must  be 
positive  if  eik2X  is  to  decay  and  produce  bound  modes  outside  the 
beam.  Also,  from  equation  (41), 

k.2  .  a2  _  k2  (47) 


so 


(48) 


Now  equation  (48)  satisfies  the  assumption  that  k,  is  real  only 

when 


L 


>  i 


(49) 


15 

i 

i 

» 


One  can  rewrite  equation  (46)  as 


tan 


b  ,  ,  ,..1/2 
-  fa2  -  k2l 
2 


ki 


Y3(k  -  aB)2 


—11/2  \ 

-'}  I 


(50) 


Y3(k  -  aB)2 


-  1 


1/2 


For  a  given  set  of  parameters  such  as  electron  density, 
velocity,  and  beam  width,  there  are  successive  values  of  a  that  satisfy 
equation  ( 50 ) .  Each  of  these  values  corresponds  to  a  mode  whose  phase 
velocity  is  a»/a.  Further  analysis  of  these  modes,  including  the 
calculation  of  these  modes  for  a  certain  set  of  conditions,  follows  in 
section  4.1. 

3.2  Cylindrical  Beam 

The  discussion  in  section  3.1  (up  to  eq  30)  applies  to  any 

restricted  beam  where  V ’  •  V  -  0.  Thus,  the  examination  of  the 

x  y 

cylindrical  beam  begins  with  equation  (30).  This  time,  in  order  to 
impose  cylindrical  symmetry  on  the  problem,  we  will  let  3/30  *  0.  In 
this  case,  the  component  forms  of  equations  (1)  and  (2)  become, 
recalling  that  F(r,0,z)  =  F(r,0)eiaz, 

ikHr  -  -iaE0  ,  -ikEr  -  iaHg  ,  (51a,b) 
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I 


iim 


Here  again,  to  get  TM  modes,  H_  is  set  to  0.  Solving  for  E,  in 

z  r 

equation  (51b)  and  substituting  the  solution  into  equation  (52a),  we 
have 


H  =  ik--  i_  E 

6  .9  2  3*  2 

k *■  -  or 


(54) 


and  equation  (30)  can  be  written  as 


flHb-i) 


E  =  0 
z 


r  <  a 


(55) 


E  =  0 
z 


r  >  a 


(56) 


where  and  k^2  are  defined  as  before  and  a  is  the  radius  of  the 


beam.  These  are  Bessel's  equations  and  have  solutions. 


M0(klr)  '  r  < 


(57) 


BKQ(k^r)  ,  r  >  a  , 


where  A  and  B  are  constants. 
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Using  equations  (51b)  and  (54)  and  the  Bessel  recursion 
formula,  **  we  find  for  r  <  a 


~~  [x  nj  (x)]  »  -X  nj  (x) 
dx  1  n  1  n+1 


•  “o(V) 


(58) 


ictk 


;r  "  ’  T T  M1^klr)  ' 


k2 

2 


(59) 


ikk 

H  - - 

6  v2 


-  AJ^^r) 


(60) 


and  we  find  for  r  >  a 


Ez  "  BK0(k2r)  ' 


(61) 


iak2 

k2 

ikk' 

H.  -  *  ~r  “iM  • 


(62) 


(63) 


Both  sets  of  solutions  must  agree  at  r  ■  a  so  that 


Abramovitz  and  I.  A.  Stegun,  Handbook  of  Mathematical  Functions, 
National  Bureau  of  Standards  (November  1970),  355-434 . 
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W0(v)  "  BK0^k2a^  * 


and  must  be  continuous,  giving 


Mc1J1(k.|a)  =  Bk^K^k^a) 


Dividing  equation  (65)  by  equation  (64)  yields 


k1  J1^kia) 


k^  ~ov~r 


Vk2aJ 


Once  more,  the  values  of  a  that  satisfy  this  relation  represent 
modes  that  can  propagate  down  the  beam.  Section  4  contains  analyses  of 
equations  (50)  and  (66),  including  the  determination  of  some  modes  for 
both  cases. 


4.  ANALYSIS  OF  BEAMS 


4.1  Analysis  of  Ribbon  Beam 


For  a  ribbon  beam  of  thickness  b,  the  modes  that  can  propagate 
are  determined  by  solving  equation  (50).  So  in  terms  of  kj  and  k^. 


M-5  • 


Let  us  begin  by  defi».  \ng  g  and  n  such  that 


k2 

9  “  — ,  n  *  1  -  r  $  •  (68a, b) 

y3k2  * 


Equation  ( SO )  becomes 


To  justify  some  of  the  following  approximations,  it  is 
appropriate  to  define  what  we  consider  realistic  values  for  the 
parameters  of  equation  (69).  The  current  state  of  electron  beam  devices 
makes  relativistic  beams  with  currents  of  several  amperes  obtainable, 
depending  on  the  size  of  the  beam.  The  values  used  here  are  designed  to 
describe  a  system  with  a  relativistic  beam,  $  =  0.7  typically,  and  a 
current  of  about  200  mA  for  a  cylindrical  beam  of  0.5-cm  radius.  This 
current  yields  an  electron  density,  nQ,  of  about  7.5  x  107 
electrons/cm3.  These  same  numbers  are  used  for  the  ribbon  beam. 
Because  of  the  interest  of  producing  power  sources  with  operating 
frequencies  between  200  and  300  GHz,  we  consider  an  output *wave length  of 
1  mm.  In  summary,  the  calculations  are  carried  out  with  the  following 
values : 


X  ■  0.1  cm 

b,  a  ■  0.5  cm 


20 


i 


n  *  7.5  x  107/cm3 
o 


3  «  0.7 

which  correspond  to  198  raA  for  the  cylindrical  beam  or  252  mA/cm2  for 
the  sheet  beam  (recall  that  we  did  not  specify  a  boundary  in  the  y- 
direction) . 


To  continue  with  the  analysis,  to  get  real  arguments  for 
equation  (69),  n  must  be  between  +g  and  -g.  For  the  values  of  k^,  k, 
and  6,  g  is  fairly  small  (about  1.5  x  10~4),  so  we  first  let 


(70) 


and  then  make  the  approximation  n  *  0  in  the  expression 


[(1  -  n)2  -  e2]l/2  ~  1/Y  • 

Mow  equation  (69)  can  be  written  as 


Finally,  we  further  simplify  the  argument  of  tan  by  letting 


Z 


(72) 
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so  that 


Z 


tan  Z  = 


kb 

2BY 


(73) 


The  roots,  Zffl,  of  equation  (73)  are  not  difficult  to  find.  A 
brief  glance  at  the  function  z  tan  Z  shows  that  there  must  be  values  of 
Z  in  every  it  interval  that  satisfy  equation  (73).  Figure  t  illustrates 
this  result  for  the  case  where  g,  b,  and  £3  have  the  values  chosen 
earlier.  Only  the  positive  portion  of  Z  tan  Z  has  been  plotted.  Table 
1  contains  the  calculated  values  of  Z  corresponding  to  the  first  11 
modes. 


The  third  column,  Z  tan  Z  -  kb/2BY»  gives  some  indication  of 
the  accuracy  of  the  roots. 

Combining  equations  (70)  and  (72)  gives 


(74) 


Thus,  for  every  value  of  Z  there  are  two  values  of  n  to  be 
considered  in  the  exact  case,  that  is,  when  n  is  no  longer  approximated 
as  zero  in  [(1  -  r\)2  ~  B2]1^2.  In  this  way,  two  values  of  a  are 
generated  corresponding  to  "fast"  and  "slow"  solutions  of  equation 
(50).  Table  2  contains  the  fast  and  slow  propagation  constants  for  the 
first  11  modes. 


The  fast  solutions,  corresponding  to  smaller  a's,  are  found  by 
taking  the  positive  case  in  equation  (74).  Slow  solutions  are  generated 
by 'using  the  negative  case. 


Figure  1 .  Behavior  of  2  tan  Z  as  function  of  Z 
demonstrating  periodicity  of  values  of  Z  that 
satisfy  equation  (73):  X  =  0.1  cm,  b  =  0.5  cm, 
and  0  =  0.7. 


TABLE  1.  ANALYSIS  OF  RIBBON  BEAM, 
APPROXIMATE  ROOTS:  X  -  0.1  cm, 
b  «  0.5  cm,  nQ  -  7.5  x  107/co3, 
AND  0  -  0.7 


Mode  number,  m 

Z 

Z  tan  Z  -  kb/28X 

0 

1.47878 

1.144  *  10~5 

1 

4.44199 

9.5367  x  10"6 

2 

7.42034 

6.485  x  10~5 

3 

10.4191 

o 

© 

4 

13.4393 

9.5367  x  10-6 

5 

16.4794 

-2.0981  x  10-S 

6 

19.5365 

2.4795  x  10”5 

7 

22.6078 

5.722  x  10-6 

8 

25.0905 

3.6239  x  1C"5 

9 

28.7824 

3.2424  x  10*5 

10 

31.8817 

2.6703  x  10~5 
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TABLE  2.  ANALYSIS  OF  RIBBON  BEAN, 

EXACT  ROOTS:  A  »  0.1  ca,  b  -  0.5  cat, 
n3  •  7.5  x  107/cm3,  AND  6  -  0.7 


Mode  number,  m 

a  fast 

a  slow 

0 

89.7459 

89.7739 

1 

89.7464 

89.7734 

2 

89.7472 

89.7726 

3 

89.7481 

89.7716 

4 

89.7491 

89.7706 

5 

89.7501 

89.7697 

6 

89.7510 

89.7688 

7 

89.7518 

89.7680 

8 

89.7525 

89.7673 

9 

89.7531 

89.7667 

10 

89.7536 

89.7662 

One  can  see  from  figure  2 
that,  as  the  mode  number  increases, 
the  difference  between  a  fast  and  a 
slow  decreases.  In  fact,  ot  fast  and 
o  slow  converge  to  the  value  of  o 
that  would  yield  a  phase  velocity 
equal  to  the  particle  velocity  of 
the  beam;  that  is,  a  =  k/8. 

It  is  possible  to 
investigate  the  ratio  of  energy 
flowing  outside  the  beam  to  that 
flowing  inside  for  various  modes  by 
integrating  the  average  Poynting 
vector  over  a  cross  section  of  the 
beam.  Thus,  we  want  to  evaluate 


f  Sav  •  da  =  ^  f  Re(E  x  H*) 
Js  Js 


da 


(75) 


Since  we  are  interested  in  energy  flowing  in  the  z-direction 


and  Hx  =  Hz  =  Ey  =  0, 


W  =  -  /  RefE  H*)  dx  dy  , 

2  Js  X  y 


H*  =  - — —  [~k  A  sin  (k  x)e  iOZ] 

y  V2  -  *2  L  1  1  ; 


E  = 


io 


[-k^A  sin  (k1x)ei0t“]  > 


X  ,2  ? 

k*  -  a* 


(76) 


(77). 


(78) 
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Outside  the  beam 


k  (ikiBe-k2xe~iaz)  , 


H*  = 

y  k2  -  a2  2 


(81) 


x  2  2 

k^  -  or 


^ -  (ik'Be“k2xeiaz)  ; 


(82) 


thus,  the  energy  flowing  outside  the  beam,  Wn,  is 


»o  -r  r 


-  v.2-2  -2k2x 


0  ~b/2  (k2  -  a2)2 


k' zBze  z  dx  dy 


(83) 


akk'B2 


2  e-k2b 


2(k2  -  a2)2 


(84) 


From  equation  (45),  we  have 


v  (hl\  icjb/2 
k;  Sln  \  2  /  6 


(85) 


so  equation  (84)  becomes 


ak2kA2 

- 1 -  sin2 

2k^(k2  -  a2}2 


(86) 
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In  equations  (80)  and  (86),  A  is  determined  by  the  manner  in  which  the 
beam  is  perturbed  from  its  ideal  steady-state  condition.  Such 
perturbations  are  nearly  inevitable,  but  very  difficult  to  predict.  We 
can  nevertheless  consider  the  ratio  of  these  energies  independent  of 
A.  Thus,  after  cancellation  of  the  common  constants. 


(87) 


(88) 


Hereafter,  W0/(WQ  +  Wj)  is  referred  to  as  Q.  Table  3  contains  values 
for  Q  for  the  first  10  modes. 

Figure  3  compares  the  energy  ratios  for  these  modes.  For  the 
first  mode,  at  least,  the  slow  waves  are  more  efficient,  coupling  a 
larger  percentage  of  the  energy  outside  of  the  beam.  Even  with  the 
first  mode,  however,  Q  is  less  than  1.7  percent.  For  both  fast  and  slow 
waves,  the  efficiency  is  monotonically  decreasing  with  increasing  mode 
number. 
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TABLE  3.  ANALYSIS  OF  RIBBON  BEAM,  EFFICIENCY: 
A  “  0.1  cm,  b  -  0.5  cm,  nQ  “  7.5  x  107/cm3, 
AND  0  -  0.7 


Mode  number,  m  Q  fast  Q  slow 


0 

1.57364 

X 

10-2 

1.69584 

X 

10”2 

1 

1.51604 

X 

1 0”2 

1.50324 

X 

10”2 

2 

1.32180 

X 

10"2 

1.32568 

X 

10"2 

3 

1.12045 

X 

10”2 

1.17251 

X 

10“2 

4 

9.31548 

X 

10"3 

9.15241 

X 

10"3 

5 

7.69732 

X 

10‘3 

7.80541 

X 

10"3 

6 

6.61363 

X 

10"3 

6.10717 

X 

10-3 

7 

5.55854 

X 

10-3 

5.67797 

X 

10-3 

8 

4.62158 

X 

10-3 

4.74057 

X 

10-3 

9 

3.46832 

X 

10-3 

3.58016 

X 

10-3 

Figure  3.  Efficiency  as  function  of 
increasing  mode  number  for  both  fast 
and  slow  modes:  X  =  0.1  cm,  b  =  0.5  cm, 
nQ  =  7.5  x  107/cm3,  and  6  =  0.7. 
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4.2  Analysis  of  Cylindrical  Beam 


For  a  cylindrical  beam  of  radius  a,  we  derived  the  relations 

k1  JlOia)  *lfk2aj 
k2  J0lkia)  *  Vk2a^  ' 


where  and  k^  are  defined  as  before. 


If  we  again  define  g  and  n  as  in  equations  (68a)  and  (68b), 

then 


As  with  the  ribbon  beam,  we  let  z  =  (g2/n2  -  l)1^2  and  make  the 
approximation  n  ~  0  in  [(1  -  n)2  -  02 ] ^2 .  Then  equation  (90)  becomes 


■)  «,(& 
1  ».(£  ■>  «.$) 


(91) 


From  the  previous  procedure  and  equation  ( 72 ) , 


Z 


ka 

0Y 


z 
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so  that  we  have 


J^Z) 

JQ(Z)  '  0Y 


ka  Ki(gY) 


KoW 


(92) 


The  evaluation  of  this  expression  is  not  too  difficult, 
especially  when  one  recognizes  that  the  right-hand  side  of  equation  (92) 
is  constant.  Now,  Jg(Z)  and  J.,(Z)  have  zeros  occurring  with  the  same 
frequency  in  Z.  Thus,  we  would  expect  J^(Z)/JQ(Z)  to  reach  any  value 
between  -**>  and  +»  once  between  each  zero  of,  say,  Jq(Z).  This  range  is 
illustrated  in  figure  4. 


Figure  4.  Behavior  of  Z[ ( Z)/Jg (Z) ]  as  function  of  Z 
showing  periodicity  of  values  of  Z  that  satisfy 
equation  (92):  X  =*  0.1  cm,  a  *  0.5  cm,  and  different 
values  of  0. 
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By  Henkel's  asymptotic  expansion**  for  large  Z, 


J^( Z)  »  4irZP(v,Z)  cos  x  -  Q(v,Z)  sin  x 


x  * 


where 


P(v,Z) 


( 4v2  -  1)(4v2  -  9) 
21 (8Z)2 


+  (4v2  -  1)(4V2  -  9)(4v2  -  25) (4V2  -  49) 
21  ( 8Z)** 


and 


Q(v,Z) 


( 4v2  -  1) 
8Z 


( 4v2  -  1 ) (4v2  -  9 ) ( 4v2  -  25)  + 
31 (8Z)2 


9 


we  can  write*  keeping  only  the  first  order  terms. 


J^Z) 

V*7 


*  tan 


Z  large 


(93) 


which  explains  the  tangent-like  behavior  in  figure  3. 


**#.  Abr amowitz  and  I.  A.  stegun,  Handbook  of  Mathematical  Functions, 
National  Bureau  of  Standards  (November  1970),  355-434. 


In  the  actual  calculation  of  equation  (92),  a  short  power- 
series  routine  was  used  for  small  values  of  Z.  Machine  limitations 
reduced  accuracy  for  Z  >  10,  so  we  resorted  to  Hankel's  expansion  for 
larger  values.  It  was  found,  however,  that  the  second  order  terms  were 
necessary  to  achieve  the  desired  accuracy  for  Z  near  10.  Table  4 
contains  some  results  for  realistic  parameters. 


TABLE  4.  ANALYSIS  OF  CYLINDRICAL  BEAM,  APPROXIMATE 
ROOT'S <  X  -  0.1  cm,  a  -  0.5  cm,  B  -  0.7,  AND 
X  •  ka/BX 


Mode  number,  m 

Zm 

ZJ, (Z)/J0(Z) 

XK^ ( X )/Kq ( X) 

0 

2.33219 

32.5468 

32.5468 

1 

5.35454 

32.5470 

32.5468 

2 

8.39744 

32.5460 

32.5468 

3 

11.4489 

32.5469 

32.5468 

4 

14.5062 

32.5476 

32.5468 

5 

17.5699 

32.5468 

32.5468 

6 

20.6397 

32.5472 

32.5468 

7 

23.7156 

32.5468 

32.5468 

8 

26.7972 

32.5468 

32.5468 

9 

29.8843 

32.5469 

32.5468 

10 

32.9763 

32.5474 

32.5468 

11 

36.0728 

32.5467 

32.5468 

As  with  the  ribbon  beam,  we  used  equation  (74)  and  plugged  back 
into  the  exact  expression,  equation  (90).  Table  5  lists  the  fast  and 
slow  roots.  These  results  are  plotted  in  figure  5.  Once  again,  the 
fast  and  slow  propagation  constants  are  converging  to  the  phase  velocity 
.of  a  wave  travelling  at  the  particle  velocity  of  the  beam.  This  time, 
however,  they  are  approaching  this  value  at  a  slower  rate  than  with  the 
ribbon  beam. 


r 


The  calculation  using  the  average  Poynting  vector  is  a  little 
more  complicated  for  this  geometry,  but  it  is  still  relatively 
straightforward.  As  before,  the  expression  to  evaluate  is  equation 
(75), 


jf  Kv  *  ■  i  j[ Re^ x 


5*) 


da 


and  again  we  are  concerned  with  energy  travelling  in  the  z-direction  and 
with  TM  modes,  so 


W 


J Re^Hj^r  dr  d0 


(94) 


For  r  <  a. 


ikk 

He  “  AJi(V) 

*2 


(95) 


E 

r 


(96) 


which,  when  substituted  into  equation  (94),  gives 


W 


I 


Jl(*1r)r  dr  * 


(97) 


Similarly,  for  r  >  a 


ikk' 


h; 


bk 


0.2  2  ^ 
k^  -  or 


M 


(101) 


iak ' 


BK. 


*2  -  «2 


(k2r) 


(102) 


which  yields 


akk 


H 


2  -B2 


2(k2  -  a2)2  Ja 


n® 

Kl(*t2r)r  Jr  d6  ’ 


(  103) 


and  from  equation  (97), 


nakk' 2  i 

«- - 1 — » r  [*5(^2*)  -  Ko(®i»)K2(ki*)]  •  1 104 ' 


(k2  -  a2)2 


Here, 


so 
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As  with  the  sheet  beam,  these  expressions  for  the  -energy  are 
not  useful  without  some  way  to  specify  A,  but  the  ratio  of  wq/(wo  +  Wj) 
is  independent  of  A.  We  have,  then. 


where  the  argument  of  Jg  and  is  k-^a  and  that  of  Kg  and  K1  is  k^a. 

Table  6  contains  some  calculated  values  for  wq/(wq  +  Wj )  (which 
is  called  Q) .  These  values  are  plotted  in  figure  6.  These  efficiencies 
are  roughly  twice  as  large  as  those  of  the  ribbon  beam  and  have  a 
smoother  behavior  as  the  mode  number  increases.  Consistently  for  any 
mode  number,  slow  waves  couple  out  more  energy. 


TABLE  6.  ANALYSIS  OF  CYLINDRICAL  BEAM,  EFFICIENCY: 
A  -  0.1  cm,  a  «  0.5  cm,  nQ  «  7.5  *  107/cm3, 

AND  8-0.7 


Mode  number,  m  Q  fast  Q  slow 


0 

3.20058 

X 

10" 2 

3.32486 

X 

10~2 

1 

3.12990 

X 

10*2 

3.25074 

X 

1  o"2 

2 

3.01073 

X 

10"2 

3.12456 

X 

1  o"2 

3 

2.85307 

X 

10"2 

2.95827 

X 

10"2 

4 

2.67009 

X 

10-2 

2.76532 

X 

10"2 

5 

2.47362 

X 

IQ-2 

2.55864 

X 

10"2 

6 

2.27388 

X 

10-2 

2.34890 

X 

1 0*2 

7 

2.07866 

X 

10"2 

2.14423 

X 

10“2 

8 

1.89319 

X 

10"2 

1.95021 

X 

10_2 

9 

1.72068 

X 

1 0~2 

1.77009 

X 

10"2 
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MODE  NUMBER,  m 


Figure  6.  Efficiency  as  function  of  mode  number  for  both  fast 
and  slow  modes:  \  =  0.1  cm,  a  =  0.5  cm,  nQ  =  7.5  x  107/cm3, 
and  g  =  0.7. 

5 .  CONCLUSIONS 

Except  in  the  ideal  case  of  a  perfectly  unperturbed  electron  beam, 
bound  electromagnetic  modes  will  be  propagating  down  a  beam  of  charged 
particles  in  free  space.  The  behavior  of  these  modes  is  quite  similar 
for  both  the  beam  geometries  considered.  In  fact,  for  both  types,  the 
dispersion  relations  produce  similar  values  for  the  propagation 
constant.  This  similarity  is  not  surprising  since  the  dispersion 
relation  for  the  cylindrical  beam  can  be  characterized  by  a  tangent 
function  for  larger  arguments  and  the  ribbon  beam  dispersion  relation  is 
in  fact  a  tangent  function. 

It  was  interesting  to  examine  the  percentage  of  energy  flowing 
outside  the  beam  because  if  an  attempt  to  use  these  modes  were  made,  one 
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would  almost  certainly  try  to  couple  to  them  outside  the  beam.  For  our 
parameters,  the  maximum  efficiency  for  a  single  mode  was  little  more 
than  3  percent.  Although  the  numbers  were  very  similar,  the  slow  modes 
couple  out  slightly  more  of  the  total  electromagnetic  energy. 

A  more  difficult  problem  is  to  be  able  to  predict  and  characterize 
perturbations  on  the  beam  in  the  laboratory.  In  addition  to  increasing 
our  understanding  of  the  beam's  behavior,  it  might  enable  us  to  excite 
the  beam  for  some  desired  result. 

The  problems  considered  here  should  not  be  taken  as  a  result,  but  as 
a  step  toward  the  solution  of  more  complicated  problems.  Of  immediate 
technological  interest  is  the  inducement  of  Cerenkov  radiation  in  the 
millimeter  wave  range.  The  regions  outside  the  beam  considered  here 
could  be  occupied  by  dielectric  or  partially  filled  with  dielectric. 
The  ribbon  beam  case  considered  here  could  be  applied  to  the  dielectric 
configuration  considered  previously.*  The  calculation  of  the  fraction 
of  the  total  energy  propagated  outside  the  electron  beam  could  be 
repeated  for  the  Cerenkov  radiation  induced  in  the  dielectric  slab. 
Such  calculations  could  be  used  to  optimize  parameters  of  the  system  for 
the  production  of  Cerenkov  radiation  of  a  particular,  predetermined 
frequency. 
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ATTN  TECHNICAL  LIBRARY 
FT  MONMOUTH,  NJ  07703 

DIRECTOR 

NIGHT  VISION  6  ELECTRO-OPTICS  LABORATORY 
ATTN  TECHNICAL  LIBRARY 
ATTN  DELNV-L ,  DR.  R.  BUSER 
FT  BELVOIR,  VA  22060 

DIRECTOR 

US  ARMY  SIGNALS  WARFARE  LABORATORY 
VINT  HILL  FARMS  STATION 
ATTN  TECHNICAL  LIBRARY 
WARRENTON,  VA  22186 

UNDER  SECRETARY  OF  DEFENSE  RE  6  ENG 
ATTN  DEP  UNDER  SECRETARY  RE  6  ADVANCED 
TECHNOLOGY 

WASHINGTON,  DC  20301 

OFFICE  OF  ASSISTANT  SECRETARY  OF  THE  ARMY 
RE.  DEV  £  ACQ 

ATTN  PRINCIPAL  DEPUTY  ASSISTANT  SECRETARY 
(RAD) 

ATTN  DIR  OF  ARMY  RES 

DAMA-ARZ-A,  DR.  M.  E.  LASSER 
DEPARTMENT  OF  THE  ARMY 
WASHINGTON,  DC  20310 

COMMANDER 

US  ARMY  MATERIEL  DEV  £  READINESS  COMMAND 
ATTN  DRCLDC.  J.  BENDER 
5001  EISENHOWER  AVE 
ALEXANDRIA,  VA  22333 

COMMANDER 

US  ARMY  RESEARCH  OFFICE  (DURHAM) 

ATTN  DR.  R.  LONTZ 
ATTN  DR.  J.  SUTTLE 
ATTN  DR.  W.  FLOOD 
PO  BOX  12211 

RESEARCH  TRIANGLE  PARK,  NC  27709 

COMMANDING  OFFICER 
US  ARMY  FOREIGN  SCIENCE  £  TECHNOLOGY 
CENTER 

FEDERAL  OFFICE  BLDG 
ATTN  TECHNICAL  LIBRARY 
220  7TH  STREET,  NE 
CHARLOTTSVILLE,  VA  22901 
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COMMANDER 

US  ARMY  MISSILE  COMMAND 
ATTN  TECHNICAL  LI8RARY 
ATTN  CHIEF,  RESEARCH  DIRECTORATE 
REDSTONE  ARSENAL,  AL  35809 

COMMANDER 

US  ARMY  ARMAMENT  RES  S  DEV  COMMAND 
ATTN  TECHNICAL  LIBRARY 
DOVER,  NJ  07801 

COMMANDER 

US  ARMY  COMMUNICATIONS  S  ELECTRONICS 
MATERIEL  READINESS  COMMAND 
ATTN  TECHNICAL  LIBRARY 
FT  MONMOUTH,  NJ  07703 

CHIEF  OF  NAVAL  RESEARCH 
OEPT  OF  THE  NAVY 
ATTN  TECHNICAL  LIBRARY 
ARLINGTON,  VA  22217 

SUPERINTENDENT 
NAVAL  POSTGRADUATE  SCHOOL 
ATTN  LIBRARY,  CODE  2124 
MONTEREY.  CA  93940 

DIRECTOR 

NAVAL  RESEARCH  LABORATORY 
ATTN  2600.  TECHNICAL  INFO  DIV 
ATTN  DR.  V.  L.  GRANATSTEIN 
ATTN  DR.  PHILIP  SPRANGLE 
WASHINGTON.  DC  20375 

COMMANDER 

NAVAL  SURFACE  WEAPONS  CENTER 
ATTN  DX-21,  LIBRARY  DIV 
DAHLGREN,  VA  22448 

COMMANDER 

NAVAL  SURFACE  WEAPONS  CENTER 
ATTN  E-40,  TECHNICAL  LIBRARY 
WHITE  OAK,  MD  20910 

COMMANDER 

NAVAL  WEAPONS  CENTER 
ATTN  TECHNICAL  LIBRARY 
CHINA  LAKE,  CA  93555 

COMMANDER 

HQ  AIR  FORCE  SYSTEMS  COMMAND 
ANDREWS  AFB 

ATTN  TECHNICAL  LIBRARY 
WASHINGTON,  OC  20334 


SECRETARY 

DEPARTMENT  OF  ENERGY 

G-231* 

ATTN  DR.  M.  JOHNSON 
ATTN  DR.  M.  MURPHY 
ATTN  DR.  J.  WILLIS 
WASHINGTON,  DC  20545 

COLUMBIA  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 

ATTN  DR.  S.  P.  SCHLESINGER 

NEW  YORK,  NY  10027 

STANFORD  UNIVERSITY 
SLAC 

ATTN  DR.  JEAN  LEBACQ2 
STANFORD,  CA  94305 

DARTMOUTH  COLLEGE 
PHYSICS  DEPARTMENT 
ATTN  DR.  JOHN  WALSH 
HANOVER,  NH  03755 

US  ARMY  ELECTRONICS  RESEARCH  6 
DEVELOPMENT  COMMAND 
ATTN  TECHNICAL  DIRECTOR.  DRDEL-CT 

HARRY  DIAMOND  LABORATORIES 

ATTN  CO/TD/TSO/D I V I S ION  DIRECTORS 

ATTN  RECORD  COPY,  81200 

ATTN  HDL  LIBRARY,  81100  (2  COPIES) 

ATTN  HDL  LIBRARY,  81100  (W00D8RIDGE) 

ATTN  TECHNICAL  REPORTS  BRANCH,  8I3OO 

ATTN  CHAIRMAN,  EDITORIAL  COMMITTEE 

ATTN  LEGAL  OFFICE.  97000 

ATTN  CHIEF,  00210 

ATTN  CHIEF,  11000 

ATTN  CHIEF,  13000 

ATTN  CHIEF,  15000 

ATTN  CHIEF,  21300 

ATTN  CHIEF,  22000 

ATTN  CHIEF,  95000 

ATTN  H.  BRANDT,  22300 

ATTN  0.  MEYER,  22800 

ATTN  Z.  G.  SZTANKAY,  13300 

ATTN  N.  KARAY I  AN  I S,  13200 

ATTN  J.  MCGARRITY,  22300 

ATTN  J.  NEMARICH,  13300 

ATTN  J.  SILVERSTEIN,  13300 

ATTN  R.  MORRISON,  13500 

ATTN  J.  SOLN,  22300 

ATTN  R.  CURNUTT,  I33OO 

ATTN  R.  LEAVITT,  13200 

ATTN  G.  D.  DOCKERY,  13300  (10  COPIES) 

ATTN  C.  MORRISON.  13200  (10  COPIES) 


